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We study transport in biased Josephson junctions made of nanowires with strong spin-orbit cou- 
pling, as it transitions into a topological superconducting phase for increasing Zeeman field B. 
Despite the absence of a fractional steady-state ac Josephson current in the topological phase, the 
dissipative Multiple Andreev Reflection current Idc at different junction transparencies is particularly 
revealing. It exhibits unique features related to topology, such as the gap inversion, the formation 
of Majorana bound states, and fermion-parity conservation. In contrast, and rather surprisingly, 
the critical current I c does not vanish at the critical point where the superconducting gap vanishes, 
exhibiting a discontinuous dI c /dB instead. These results demonstrate the feasibility to probe the 
formation of Majorana states without the need of phase sensitive or noise-related measurements. 



Semiconducting nanowires with a strong spin-orbit 
coupling in the proximity of s-wave superconductors and 
in the presence of an external magnetic field are a promis- 
ing platform to study Majorana physics. Theory predicts 
that above a critical Zeeman field B c = y/ fi 2 + A 2 , de- 
fined in terms of the Fermi energy /i and the induced 
s-wave pairing A [TJ |2], the wire undergoes a topologi- 
cal transition into a phase hosting zero energy Majorana 
bound states (MBS) at the ends of the wire. Recent 
experiments have reported conductance measurements 
that support the existence of such MBSs at normal- 
superconductor (NS) junctions in InSb [5J 0] and InAs 
[5] semiconducting nanowires. The main result of these 
experiments is an emergent zero-bias anomaly (ZBA) in 
differential conductance dl/dV measurements as an ex- 
ternal Zeeman field B increases. In this context, the ZBA 
results from tunnelling into the MBS [MS]- Although 
these experiments are partially consistent with the Ma- 
jorana interpretation [5], other mechanisms that give rise 
to ZBAs, such as disorder [I0HI2] and Kondo physics [15] . 
cannot be completely ruled out. Moreover, the expected 
superconducting gap inversion was not observed. 

Stronger evidence could be provided by the di- 
rect observation of non-Abelian Majorana interfer- 
ence (braiding) |14) . or by transport in phase-sensitive 
superconductor-normal-superconductor (SNS) junctions. 
The latter approach, which typically involves the mea- 
surement of an anomalous "fractional" 47r-periodic ac 
Josephson effect [T5HT7] is in principle much less de- 
manding than an interference experiment. Realistically, 
however, such fractional current may be difficult to de- 
tect, since dissipation is expected to destroy it in the 
steady state regime. Although it has been shown that 
the anomalous Josephson current survives in the dynam- 
ics, such as noise and transients [TBl420| . simpler experi- 
mental probes into MBS are extremely desirable |21) . 

Here we propose the multiple Andreev reflection 
(MAR) current in voltage-biased SNS junctions [2"21 125] . 
Idciy), as an alternative, remarkably powerful, yet sim- 
ple tool to study the topological transition of junctions 
made of semiconducting nanowires. This is made pos- 



sible by the direct effect that gap inversion, MBS for- 
mation and fermion-parity conservation have on at 
various junction transparencies. For vanishing normal 
transparency T/v, Idc{V) traces the closing and reopen- 
ing of the superconducting gap at B c , A off ~ \B — B c \. 
This gap inversion can be shown to be a true topological 
transition by tuning the junction to perfect transparency 
T N = 1. The dissipative MAR current I dc {V -> 0) in 
this limit shows parity conservation effects in the pres- 
ence of MBS that are just as significant as a fractional ac 
Josephson effect would be, but do survive in the steady 
state limit. Moreover, the detailed dependence of MAR 
as a function of Tjv has the fundamental advantage over 
transport spectroscopy in NS junctions in that it contains 
information about the peculiar dependence of MBS hy- 
bridization with superconductor phase difference (j) and 
Tjv, despite not requiring any external control of </>. Sim- 
ilarly, we show that another important phase-insensitive 
quantity, the critical current I c , has an anomalous be- 
haviour across the topological transition. Despite the 
expectation that it should vanish with the supercondut- 
ing gap at B Cl we show that it remains finite due to a 
significant continuum contribution, revealing the gap in- 
version only as discontinuity in dI c /dB. 

Nanowire spectroscopy. — A single one-dimensional 
(ID) nanowire in the normal state is modelled in terms of 
its effective mass m*, spin-orbit coupling a so and Fermi 
energy /z, Zeeman splitting B = gfigB, where g is the 
semiconductor's g-factor and B is the magnetic field along 
the wire. The wire is close to depletion, so \x is not very 
large. It's normal Hamiltonian reads (<Xj are spin Pauli 
matrices) 

P 2 

H o = K~^r _ M + Bo- x + a so o- y p. 
2m* 

In the presence of an induced s-wave superconducting 
pairing A (assumed real), the system is modeled by the 
Nambu Hamiltonian 
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FIG. 1. (Color online) (a) Lowest bands of a B = 0.5B C 
nanowire, with (dashed) and without (solid) pairing A. (b) 
Evolution of bands with Zeeman field B. The gap A_ closes 
at B — B c , while A+ does not. (c) Local density of states at 
the junction for perfect normal transparency Tn = 1, which is 
peaked at the energy e±(4>) of Andreev (quasi)bound states. 
The minimum 8^ of e.±(<t>) for B < B c falls to zero at perfect 
transparency as /i/A is increased (Andreev approximation). 
Panel (d) shows the Tn = 1, B — <y = 5^. It falls as 
Sir = ciA 2 /([i + AC2) for some ci,2 > 0. 



The essential ingredient for a topological supercon- 
ductor is an effective p-wave pairing acting on a single 
("spinless") fermionic species, the so-called Kitaev model 
|15j . Spin-orbit coupling splits nanowire states into two 
subbands of opposite helicity at B = 0. At finite B, 
these two subbands, which we label 4- and — [black and 
gray lines in Fig. [IJa)], have spins canted away from 
the spin-orbit axis. The s-wave pairing A, expressed 
in the ± basis, takes the form of an intraband p-wave 
Ap +,/ (p) = ±ipAa so /^B 2 + (a so p) 2 , plus an inter- 
band s-wave pairing A+~(p) = AB/^B 2 + (a BO p) 2 [24] . 
Without the latter, the problem decouples into two in- 
dependent p-wave superconductors, while A]j acts as a 
weak coupling between them. Each quasi-independent 
sector has a different gap. These two (-B-dependent) 
gaps, which we call A_ (at small momentum) and A + 
(large momentum), are shown in Fig. [lja,b). While A + 
remains roughly constant with B (for strong spin-orbit 
coupling), A_ vanishes linearly as B approaches the crit- 
ical field, A_ ss \B — B C \ 25J. This closing and reopening 
(gap inversion) signals a topological transition, induced 
by the effective removal of the — sector away from the 
low-energy problem. Below B c the nanowire is composed 
of two (weakly coupled) spinless p-wave superconducting 
species, and is therefore topologically trivial. Above B c , 
A_ is no longer a p-wave gap, but rather a normal (Zee- 
man) spectral gap already present in the normal state, 
transforming the wire into single-species p-wave super- 
conductor with non-trivial topology. This phase is host 



to pairs of MBSs at the wire ends, which are protected by 
the effective gap A c ff = Min(A + , A_). The same physics 
also emerges in the problem of superconducting helical 
edge states in spin-Hall insulators, which is spinless to 
start with [TTlfFg] . 

In a short Josephson junction made of two nanowires 
across a contact with normal transparency Tjy, an An- 
dreev bound state (ABS) is formed in each of the two 
sectors above for B < B c , while only one, associated to 
A + is formed for B > B c . Their energies e± depend on 
the superconducting phase difference 4> across the junc- 
tion, as shown in Fig. [TJc) (case Tn = 1). Due to the 
residual A s coupling between the two sectors, the A + An- 
dreev state is actually quasibound in the energy window 
A_ < e < A + . Moreover, in the non-topological phase 
(left panel) their energy does not reach zero at 4> — n, un- 
like predicted by the standard theory for Tn = 1 within 
the Andreev approximation > A |26| . The energy min- 
imum 5n does indeed vanish as /u/A grows, see Fig. [l^d). 
In the topological phase (right panel), the surviving ABS 
associated to A + arises as the hybridisation of two MBS 
across the junction. Global fermion-parity conservation 
protects the <fi = ir level crossing, giving exactly 6^ = 0. 

ac Josepshon in nanowire junctions. — Here we are con- 
cerned with the effects of the topological transition in 
the transport across a short junction of two semi-infinite 
wires. We wish to compute, for different normal trans- 
parencies Tn, the MAR current Idc{V) as a function of 
bias V across the junction, and the critical current I c 
supported at V = 0. For computation purposes, we con- 
sider a discretisation of the continuum model Eq. ([IJ for 
each wire into a tight-binding ID lattice with a small 
lattice spacing a [used also to compute Fig. [ljc,d)]. 
This transforms terms containing the momentum opera- 
tor p into nearest-neighbour hopping matrices v. Namely 



Ho = T, i 4hc i +Y: 



-h.c., where t = h 2 /2r 



and the spin structure of h and v is 



2t- n B 
B 2t- fjt 



-t 



-t 



The coupling across the junction is modelled by a hop- 
ping matrix vo = vv between the end sites of each wire, 
labeled I and r, that is suppressed by a dimensionless 
factor v € [0,1], which controls the junction's normal 
transparency at B — 0, Tn(v). 

Under a constant voltage bias V, the pairings A to 
the left and right of the junction acquire an opposite 
and time-dependent phase difference, <p(t) = 2eVt/h. 
This induces Landau-Zener transitions between the ABS 
and into the continuum, thereby developing a time de- 
pendent Josephson current with both I^c and I ac com- 
ponents. Such is the point of view exploited e.g. in 
Refs. [in] and [57]. Alternatively, the time depen- 
dence of 4>(t) can be gauged away into the hopping 
term across the junction, which then becomes Vo(t) = 



-iu!Qtr z 



E 



,i C r(T V a 



ic\ a i + h.c, where t z is the z-Pauli 



matrix in Nambu space. Driving frequency is ujq = eV/h. 
By employing Keldysh-Floquet theory [35] [55] , see Ap- 
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pendix, we obtain the stationary-state time-dependent 
current which contains all the harmonics of the ac Joseph- 
son effect, i.e. I(t) = ^3 e mUot I n . Here, we concentrate 
on the dc-current Idc = la as a function of V and T/y. 
The results for Idc{V) at small, intermediate and full 
transparency are summarised in Fig. [2^a-c) for increas- 
ing values of B spanning the topological transition (blue 
is the non-topological phase, red the topological one host- 
ing MBS). 

Tunneling regime. — As the junction's normal trans- 
parency Tjv approaches zero in the non-topological 
regime, dc-transport vanishes below an abrupt thresh- 
olds voltage V t = 2A off /e = 2A_/e (Fig. ga), blue 
curves), where, recall, A e g- = Min(A_,A + ). This well 
known result follows from the fact that there are no 
quasiparticle excitations in the decoupled wires for en- 
ergy e £ (— A e g , A e ff) if B < B c . Indeed, to second order 
in perturbation theory in z/, the MAR current takes the 
form of a convolution between A (u>) and A (u> ± eV/h), 
where Aq is the decoupled (v ~ 0) spectral density at 
each side of the junction. [The trace of Ao(e), propor- 
tional to the local density of states (LDOS), is shown in 
Fig. [3]^a)]. Hence, as B increases, the tunnelling current 
threshold follows the closing of the gap in the LDOS, 
until Vt vanishes and Idc becomes linear in small V at 
B c (black curve). As B > B c , the gap reopens, but the 
threshold is now halved to V t — A e ff/e (red curves) [50] . 
The change, easily detectable as a halving of the slope 
of the threshold dVt/dB across B c , is due to the emer- 
gence of an intra-gap zero-energy MBS in the topological 
phase [see zero energy peak in Fig. [3ja)], which opens a 
tunnelling transport channel from or into the new zero 
energy state. Moreover, as B is increased above a certain 



(a) T N =0.05 



T N =0.5 



surpasses A + , and A c g saturates at A^ 
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field B c 2) , A_ 

This is directly visible in Vt(B) as a kink at B c 
dashed and dotted lines in Fig. |2]ja)]. [51] . 

Intermediate transparency regime. — In higher trans- 
parency non-topological junctions, subharmonic MAR 
steps develop at voltages Vt/n = 2A c s/en (n = 
2,3,4,...), see Fig. [2jb). The specific profile of each 
step with V still contains information on the LDOS of the 
junction at energies around A e g. At B = 0, the power- 
law LDOS for |e| > A results in a staircase-like curve 
Idc{V) [blue line in Fig. |2jd)]. This shape is roughly pre- 
served up to B = B c . For B > B c the MAR subharmonic 
profile changes qualitatively, however. The subharmonic 
threshold voltages Vt/n are halved (since Vt is halved), 
and the MAR current profile becomes oscillatory instead 
of step-like. A blowup of the oscillations is presented in 
Fig. J2^d) (red curve), together with guidelines for the 
corresponding Vt/n in gray. 

The emergence of oscillatory MAR steps is connected 
to the formation of the zero energy MBS, specifically to 
the sharp zero energy peak in the LDOS that is a con- 
sequence of the localized nature of the MBS. This effect 
is well known from the studies of the MAR current in 
Josephson junctions containing a resonant level [52TI34) . 
such as for instance a quantum dot between two super- 






O I (V ""^ B 
-- 2A_G„\ 


SB 


2) 


- 2A + G \ 







1 

B/B c 



FIG. 2. (Color online) Time-averaged Josephson current Idc 
as a function of bias V for increasing Zeeman field B. Curves 
are offset by a constant 2A<5o / e, with Go — e 2 /h. Blue and red 
curves correspond to the non-topological (B < B c ) and topo- 
logical (B > B c ) phases respectively. Panels (a) to (c) show 
the cases of tunnelling, intermediate and full transparency. 
Panel (d) is a blowup of the low bias MAR subharmonics at 
intermediate transparency. Panel (e) shows the asymptotic 
Idc(V — > 0) at full transparency (circles), along with the de- 
pendence of the quantities 2A_(Jo and 2A+C?o with B across 
the topological transition [dashed/dotted lines, evolution also 
shown in panel (a)]. 



conductors. In the latter case, however, the oscillatory 
MAR subharmonics arise at odd fractions of 2A c fj, i.e. 
at voltages 2A c ff/(2n + I)e, instead of the A c s/en of the 
Majorana case. This difference is ultimately due to the 
fact that a quantum dot resonant level spatially localised 
within the junction cannot carry current directly into the 
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reservoirs of the external circuit, while a zero energy MBS 
(essentially half a non local fermion) can. This same sit- 
uation arises in d-wave Josephson junctions [35], which 
also exhibit oscillatory A c g/en MAR subharmonics ow- 
ing to the presence of mid gap states. 

Transparent limit. — In the limit T/v — > 1, ABS energies 
£±{<j>) [Fig. [l|c)] touch the continuum at <f) — 0. From 
the Landau-Zener point of view of the ac Josephson ef- 
fect [27], the time dependence of (f>(t) = 2eVt/h for an 
arbitrarily small V will induce the escape of any quasi- 
particle occupying an ABS into the continuum after a 
single <j)(t) cycle. A given ABS becomes occupied with 
high probability in each cycle around <j) = tt if the rate 
hd(f>(t)/dt = 2eV exceeds its energy minimum e(7r) = 8 n . 
At voltages exceeding such minimum, one quasiparticle 
is injected into the continuum per cycle, and a finite 
Idc{V > (>7r/e) arises. Below such voltage, however, the 
ABS remains empty, so that if <5j is finite, as is the case 
of a realistic non-topological junction (see Fig. [ljc)), one 
obtains Idc(V —> 0) = (valid for any transparency at 
B < B c ). This is in contrast to the conventional B = 0, 
T N = 1 result I(V -> 0) = 4A£ /e, predicted within the 
Andreev approximation (Qo = e 2 /h). 

After the topological transition this picture changes 
dramatically. The two MBS at each side of the junction 
hybridise for a given <f> into a single ABS. This seem- 
ingly innocent change has a notable consequence. Since 
fermion parity in the superconducting wires is globally 
preserved, an anticrossing at = tt, which would repre- 
sent a mixing of a state with one and zero fermions in 
the lone ABS, is forbidden. Parity conservation there- 
fore imposes 6„ = in the presence of MBS, irrespec- 
tive of Tjv or fi/A |36] . This is a true topologically pro- 
tected property of the junction, and gives rise to a finite 
Idc(V — > 0) = 2A c ffC/0i i- e - half the value expected for 
the non-topological junction in the Andreev approxima- 
tion. This abrupt change is shown in Fig. [2Jc,e). The 
Idc(V — > 0) MAR current in the transparent junction 
limit, therefore, directly probes the emergence of parity 
protected ABS crossings, formed by the fusion of the two 
Majorana modes in the junction. 

Critical current. — For finite impedance of the exter- 
nal circuit, a. V = supercurrent peak should arise su- 
perimposed onto the computed Idc [221 123 [37] . In the 
transparent limit, this peak may hinder the experimen- 
tal identification of the Idc(V 0) limit, but itself holds 
valuable information about the transition. The critical 
current I c may be computed in general by maximizing 
the V = (time-independent) current respect to <f) 
(including the contribution from the continuum). For a 
short transparent junction at B = 0, I c is maximum, and 
equal to Z° = AeA/h. Fig. ^h) shows I c for increasing 
values of the Zeeman field. Note that, in spite of the 
naive expectation that a nanowire without a supercon- 
ducting gap should not carry a finite supercurrent, this 
is not the case. At B = B Cl I c is finite, while A c ff — 
[the junction LDOS at criticality is also gapless, see Fig. 
[3]ja)]. This gapless supercurrent comes from the £+(</>) 
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FIG. 3. (Color online) (a) Local density of states at the end of 
a single nanowire in the non-topological (top), critical (mid- 
dle) and topological phase (bottom). A zero-energy Majorana 
peak appears in the latter case, (b) The critical current I C (B) 
across the the topological transition in units of 1° — AeA/h. 
The dotted line corresponds to |(A + + A_)/A for B < B c , 
and §A+/A for B > B c . 

quasi-bound Andreev state in the continuum, which con- 
tributes to the current almost as if it were a subgap ABS. 
It is thus a reasonable approximation to write I c as the 
sum of the critical current from each ABS. For B < B c , 
I c sa \lc{A + + A_)/A. The A_ contribution, however, 
should not be included for B > B c , leading to a dis- 
continuous derivative dI c /dB at B c as a signature of the 
topological transition. This simple model gives a qualita- 
tive fit to the exact numerical results [dotted line in Fig. 
[3|h)], with deviations coming from corrections to the An- 
dreev approximation, and contributions above A + . 

In conclusion, we have shown that the steady state dc- 
current in biased Josephson junctions is a flexible exper- 
imental probe into the various aspects of the topological 
superconducting transition in semiconducting nanowires. 
Tuning the junction transparency one may obtain ev- 
idence of MBS formation as conclusive as a fractional 
Josephson effect, without requiring control of the junc- 
tion phase. Moreover we have found that the critical cur- 
rent in the wire does not vanish at the transition due to 
above-gap contributions, although its derivative with B 
exhibits a discontinuity as a result of the disappearance 
of one Andreev bound state. 
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Appendix: Keldysh-Floquet theory 

We present a summary of the Keldysh-Floquet the- 
ory employed in the computation of steady-state, 
time-averaged Josephson current Id c through a biased 
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superconducting-insulating-superconducting juntion, as 
a function of bias V and junction transparency. 

Consider a mesoscopic system composed of two semi- 
infinite leads (labeled L and R), each in thermal equi- 
librium at the same temperature T and with the same 
chemical potential fi — 0. Each lead has a finite s-wave 
superconducting pairing A a , where a = L, R. A central 
system (a = S) , which may or may not be superconduct- 
ing, is coupled to both leads through operator v. In its 
Nambu form, the Hamiltonian of the system reads 



where the Nambu Hamiltonian matrix takes the general 
form 
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f h L 
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Here h a is the normal Hamiltonian for each section of the 
system. The blocks delimited by lines denote the Nambu 
particle, hole and pairing sectors. 

If we apply a left-right voltage bias V through the junc- 
tion, the BCS pairing of the leads will become time de- 
pendent, A L/R -> e ±lVt A L/R , while h L/R -> h L/R ±V/2 
(we take e = H = 1). Both these changes can be gauged 
away from the leads and into the system by properly re- 
defining cf — » cf(t) = e ±%vt ^cf. This transformation is 
done also inside the system S, thereby effectively divid- 
ing it into two, the portion with an e lVt ^ 2 phase (denoted 
Sl), and the portion with the opposite phase (denoted 
S R ). This restores H to its unbiased form, save for a new 
time dependence in hs — ► hs(Vt), which is constrained 
to the coupling between the Sl and S R , 



h s (Vt) 



hs L 



,-iVt„ } + 



n 



e* vt v Q 



It is important to note that H(t) is periodic, with an- 
gular frequency uj — V. In the steady state limit (at 
long times t after switching on the potential V) all re- 
sponse functions and observables will exhibit the same 
time periodicity (all transient effects are assumed to be 
completely damped away). In particular, the steady state 
current I(t) = I(t + 2tt/ujq), so that 



'(*) = £• 



In 



for some harmonic amplitudes /„, in general complex, 
that satisfy /„ = /*„ since I(t) is real. 

This current can be computed using the Keldysh 
Green's function formalism. |38j The standard expres- 
sion for I{t) is computed starting from the definition of 



I(t) = dtNL, where N is the total number of fermions 
in the left lead. By using Heisenberg equation and the 
Keldysh-Dyson equation, one arrives at 



I(t) = Re[J(t)} 



where 



^) = f / dt'Tv{[G r (t,t')^<(t\t) + G < (t,t')^l(t' 1 t)]r z } 

The z-Pauli matrix t z above acts on the Nambu particle- 
hole sector, 




The self energy from the left lead is defined as 
££<(t',i) = vg a L ' < (t',t)v+, where g L (t',t) = g L (t' - t) 
stands for the left lead's propagator, when decoupled 
from the system (this propagator depends only on the 
time difference since the decoupled lead is time indepen- 
dent in this gauge). We define the Fourier transform of 
g as 



dte ibJt g{t) 



The retarded propagator in Fourier space is 

s£M = r — - 

while the advanced g£(w) = [g1(cu)] + . One can com- 
pute g<{uj) = if{uj)A L (uj), where f(u) = l/(e"/ feBT + 1) 
is the Fermi distribution in the leads, and Al{oS) = 
i(g L (u>) — gKu)) is the Nambu spectral function. The 
9l/ R [ an d hi particular Al/ r (uj)\ is assumed known, or 
at least easily obtainable from h L / R and v. Finally, the 
Green functions G r (t',t) and G < (t',t) correspond to the 
propagator for the full system, including the coupling to 
the leads. (Note that, in practice, since G is inside a trace 
in J(t), only matrix elements of G inside the S portion 
of the full system are needed). The retarded G r satisfies 
the equation of motion 

[idf -H(t')]G(t',t) = 5(t'-t) 

while G < (when projected onto the finite-dimensional 
system S) satisfies the Keldysh relation 



G<(t',t) 



J dt 1 dt 2 G r (t' ,ti 
t 2 ) + 



t 2 )] G a (t 2 ,t) 



Since hs in H is time dependent, G propagators depend 
on two times; unlike Hl/r or gL/R they are not Fourier 
diagonal. Instead, we can exploit the steady-state condi- 
tion, which reads 

G(t',t) = G(t' + —,t + —) 
u> u> a 



6 



to expand the system G as a Fourier transform in i! — t 
and a Fourier series in t. We define 

oe e _ le(t -t) Gn ( e ) 
-oo 2 ^ 

The natural question is how the equation of motion is 
expressed in terms of the harmonics G n (e). It takes the 
most convenient form if we redefine G n (e) (where e is 
unbounded) in terms of the quasienergy e £ [0,huj o ], he. 
e = e + rnwo 

G m „(e) = G m _ n (e + mw ) 

This has the advantage that the equation of motion trans- 
lates to a matrix equation analogous to that of a static 
system in Fourier space 

^(e + n'wo - H n , m )G T mn {e) = S n > n 

m 

where 

H n , n = J dte^ n '- n ^H(t). 

This is known as the Floquet description of the steady 
state dynamics in terms of sidebands, which appear for- 
mally as a new quantum number n. Time dependent 
portions of H (t) act as a coupling between different side- 
bands. The effective Hamiltonian for the n-th sideband 
is the static portion of H (t), shifted by — nuo- One there- 
fore sometimes defines the Floquet "Hamiltonian" of the 



system as 

where, as before, hs n > n = J dte^ n ~ n ^ t hs(i). Likewise, 
one may define the retarded, lesser and advanced Floquet 
self-energies (a = r, <, a) as 

E£ nm (3=*»n.E2 /H (e + nwd) 

(since the leads are static, the are sideband- 

diagonal). 

The Floquet equation of motion for G^ m (e) can be 
solved like in the case of a static system. Within the S 
portion of the system, we have 

G r (e) = [e-h s --Zi(e)--£ r R (e)}- 1 

Boldface denotes the sideband structure implicit in all 
the above matrices. Similarly, the Keldysh relation takes 
the simple form 

G<(e) = G r (i) [£<(?) + S<(e)] G a (e) 

Finally, the time averaged current Idc = lo takes the 
form 

I dc = 2G / deReTi- { [G r (e)X<(e) + G<(e)EJ(e)] t z ) 

where the trace includes the sideband index. In a practi- 
cal computation, the number of sidebands that must be 
employed is finite, and depends on the applied voltage 
bias V (the tipical number scales as n max ~ vq/V). We 
employ an adaptive scheme that increases the number of 
sidebands recursively until convergence for each value of 
V. 
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